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ABSTRACT 

The Planck collaboration has recently published precise and resolved measurements of the Sunyaev- 
Zel'dovich effect in Abell 1656 (the Coma cluster of galaxies), so directly gauging the electron pressure 
profile in the intracluster plasma. On the other hand, such a quantity may be also derived from 
combining the density and temperature provided by X-ray observations of the thermal bremsstrahlung 
radiation emitted by the plasma. We find a mo del- independent tension between the SZ and the X-ray 
pressure, with the SZ one being definitely lower by 15 — 20%. We propose that such a challenging 
tension can be resolved in terms of an additional, non-thermal support to the gravitational equilibrium 
of the intracluster plasma. This can be straightforwardly included in our Supermodel, so as to fit in 
detail the Planck SZ profile while being consistent with the X-ray observables. Possible origins of 
the nonthermal component include cosmic-ray protons, ongoing turbulence, and relativistic electrons; 
given the existing observational constraints on the first two options, here we focus on the third. 
For this to be effective, we find that the electron population must include not only an energetic tail 
accelerated to 7 > 10^ responsible for the Coma radiohalo, but also many more, lower energy electrons. 
The electron acceleration is to be started by merging events similar to those which provided the very 
high central entropy of the thermal intracluster plasma in Coma. 

Subject headings: cosmic background radiation — galaxies: clusters: individual (Abell 1656) — X-rays: 
galaxies: clusters 



1. INTRODUCTION 

We are motivated by the recent, spatially resolved mea- 
surements with the Planck satellite by Ade et al. (2012) 
of the Sunyaev-Zel'dovich (1980; SZ) effect in Abell 
1656, the very rich, nearby cluster in Coma Berenices 
at z = 0.023. 

The thermal SZ effect describes how the temperature of 
crossing CMB photons is modulated by the Compton up- 
scattering off the hot electrons in the intracluster plasma 
(ICP). Its strength is given by the Comptonization pa- 
rameter y = {aT/iTieC^) J d£pe(r) integrated along line- 
of-sights across the cluster. It directly probes the elec- 
tron thermal pressure pe ~ P (2 -I- 2A)/(3 -f 5X) > 0.5 p, 
here written in terms of the ICP pressure p; with the 
cosmic hydrogen abundance X « 0.76, their ratio reads 
Pe/p — 0.52. Compared with previous observations in- 
cluding WMAP^s (see Komatsu et al. 2011, and refer- 
ences therein), the resolved Planck data improve the SZ 
probing of the cluster core and extend it into the out- 
skirts, providing a handle to the complex astrophysical 
processes in the ICP to be discussed here. 

On the other hand, the ICP pressure p « nksT/^, 
(with the mean molecular weight /i ~ 0.60) may be 
also derived from combining the density n and temper- 
ature T provided by X-ray observations of the thermal 
bremsstrahlung radiation emitted by the plasma. The 
overall trend emerging from the Planck data is toward 
a deficit in the SZ relative to the X-ray pressure, see 
Ade et al. (2012). In detail, to fit the SZ data these 
authors based on empirical formulae suggested by nu- 
merical simulations (see Nagai et al. 2007) and by X-ray 
analyses (see Arnaud et al. 2010). These formulae pro- 



vide a 'universal' pressure profile for the whole cluster 
population, or a specific version for unrelaxed clusters. 
However, when applied to fit the precise Planck SZ data 
these formulae perform inadequately, as discussed by Ade 
et al. (2012); specifically, the first version turns out to 
overshoot the data in the core, and both to appreciably 
undershoot them in the outskirts, well beyond the quoted 
uncertainties. Aimed modifications of the parameter val- 
ues in the fitting formulae, that include suppression of 
unphysical central divergencies, can improve the SZ fits 
at the cost of inconsistencies with the X-ray pressure. As 
discussed by the above authors, this is also the case with 
multiparametric fitting formulae of the type proposed by 
Vikhlinin et al. (2006) for the X-ray observables. 

Thus we are stimulated to use an orthogonal approach, 
provided by the Supermodel (SM; Cavaliere et al. 2009); 
this yields a direct link between the X-ray and the SZ ob- 
servables (see Lapi et al. 2012). We specify below how 
the SM is based on a physically motivated run of the 
entropy k{r) = kBT{r)/n'^/^{r) that underlies the ICP 
support in the gravitational potential well provided by 
the cluster dark matter, to approach hydrostatic equilib- 
rium. 

The SM improves in a number of respects upon the 
classic isothermal /3— model (Cavaliere & Fusco-Femiano 
1976), adopted by Mohr et al. (1999) and Churazov et 
al. (2012) to closely describe the central X-ray bright- 
ness profile in Coma. In fact, the SM incorporates up- 
dated, weakly cusped distributions of the dark matter 
(see Lapi & Cavaliere 2009). It also accurately describes 
the central conditions both in cool, and in non-cool core 
clusters as Coma (see Molendi & Pizzolato 2001). Fi- 
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Figure 1. Top panel: Profile of the X-ray surface brightness in the Coma Cluster; green squares refer to the XMM-Newton data by 
Churazov et al. (2012), dotted line shows their /3— model fit, while solid line illustrates our SM outcome. Bottom panel: Projected profile 
of the emission- weighted temperature in the Coma Cluster; green squares refer to the XMM-JVewton data by Snowden et al. (2008), orange 
circles to the Suzaku data by Wik et al. (2009), solid line is our SM outcome (obtained on fitting the XMM-JVewton data), with the dashed 
line representing its extrapolation into the outskirts out to the virial radius _R = 2.2 Mpc (see Churazov et al. 2012). In both panels the 
shaded areas show the associated 2 — a uncertainties. 



nally, it also describes the region beyond a few 10^ kpc 
where T dechnes outwards while n drops, so that the en- 
tropy gradually rises as shown by many X-ray data (e.g., 
Snowden et al. 2008; Cavagnolo et al. 2009; Pratt et 
al. 2010; Walker et al. 2012), and as expected on basic 
astrophysical grounds. 

2. THE SUPERMODEL VIEW 

In fact, the spherically-averaged entropy profile is to 
rise from a central level kc into an outer ramp with slope 
a toward the virial boundary R, following the pattern 
(Voit 2005; Lapi et al. 2005) 



fc(r) =k, + kn {r/Rf 



(1) 



Specifically, a central baseline kc ^ 10^ keV cm^ is pro- 
duced during the early collapse and virialization of the 
cluster core; then the intergalactic gas is condensed to 
levels n ^ 10^"^ cm^'^ in step with the general over- 
densities around 200 over the average background, while 
it is heated up to temperatures /cbT k, GM/IOR ~ 
a few keVs. These conditions imply thermal pressures 
p « 2nkBT, of order a few 10~^^ erg cm""^. The base- 
line kc may be subsequently lowered by radiative cool- 
ing, or enhanced by AGN feedback (see Lapi et al. 2003; 
Fabian 2012) and deep, energetic mergers (see Marke- 
vitch & Vikhlinin 2007; McCarthy et al. 2007). These 
processes, respectively, produce current conditions of the 
cool core kind with kc ~ 10 keV cm^, or of the non-cool 
core kind with kc > 10^ keV cm^ like in Coma (Cavaliere 
et al. 2011). 

At the outer end, the slope awl and the boundary 
value kji ~ several 10^ keV cm^ are originated as en- 
tropy is continuously produced by strong virial accretion 
shocks (for observational evidence in Coma, see Brown & 
Rudnick 2011; Markevitch 2012), and then is conserved 
and stratified as the infalling gas is compressed into the 
gravitational potential well (Tozzi & Norman 2001; Cav- 
aliere et al. 2011). 

The SM formalism simply consists in inserting the en- 



tropy run of Eq. (1) in the differential equation for hydro- 
static equilibrium of the ICP (see Cavaliere et al. 2009 
for details) ; this is easily integrated to obtain linked runs 
of the thermal pressure 
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of the density n{r) oc [p{r)/k{r)]^^^, and of the temper- 
ature T{r) oc p^/^(r) k'^^^{r). Here is the value at the 
virial boundary, while M{< r) is the gravitational mass 
distribution mainly contributed by the dark matter (see 
Lapi & Cavaliere 2009). 

In Fig. 1 we show how the SM fits the projected profiles 
of X-ray brightness Sx oc T^/^ found by Churazov et 
al. (2012), and of the emission-weighted temperature 
T measured by Snowden et al. (2008) and Wik et al. 
(2009) in Coma; the reader is referred to Fusco-Femiano 
et al. (2009) for details. From these fits we extract values 
(with their 1 — a uncertainty, and consistent with the last 
reference) of the parameters kc « 535 ±180 keV cm^, 
a « 1.3 ± 0.2, and kn w 5050 ± 225 keV cm^ specifying 
the entropy pattern in Eq. (1). 

Based on such values, Coma turns out to be an extreme 
HE (high entropy) cluster both in the core and in the 
outskirts, according to the classification used in Cavaliere 
et al. (2011). As discussed there, such conditions call for 
impacts of several energetic mergers down to the center 
that deposit energies of order lO^'* ergs, and for strong 
accretion shocks standing at the virial boundary with 
Mach numbers 10 that produce outer temperatures 

^bTr 5 keV. Both processes are in tune with the rich, 
supercluster environment that surrounds Coma. 

From Eq. (2) we find the radial pressure profile p{r) 
with no recourse to delicate deprojections. Thence we 
obtain the thermal electron pressure pc, and compute 
the profile of the corresponding SZ Comptonization pa- 
rameter y (cf. § 1)- In Fig. 2 we express our re- 




Figure 2. Profile of SZ effect toward the Coma cluster. Green squares refer to the Planck data by Ade et al. (2012), and orange circles 
to the WMAP data by Komatsu et al. (2011). Dashed line illustrates the SM outcome (smoothed on the Planck resolution scale) based on 
the fit to the X-ray data from XMM-Newton (see § 2), with the heavy shaded area representing the associated 2 — a uncertainty; dotted 
line and light shaded area illustrate the same when basing on the fit to the X-ray brightness from ROSAT data. The solid line is our 
outcome when a non-thermal contribution Sp/p 20% (or 15%) to the pressure is included in the SM (see Eqs. 2 and 3). 



suit in terms of the equivalent Rayleigh- Jeans decrement 
AT = -2 yTcMB of the CMB temperature Tcmb ~ 2.73 
K, to compare with the Planck measurements as pre- 
sented by Ade et al. (2012), see their Fig. 4. 

3. COMPARING THE SZ AND X-RAY VIEWS OF COMA 

Fig. 2 highlights a deficit in the values of |AT| as mea- 
sured by Planck, relative to those expected from the X- 
ray pressure. The discrepancy appears to be remarkably 
sharp at the center, well beyond the uncertainties budget 
presented by Ade et al. (2012). 

Such a SZ vs. X-ray mismatch goes also beyond the 
uncertainties affecting the entropy parameters from the 
X-ray fits (cf. shaded areas in Figs. 1 and 2), with effects 
damped out by the weak dependence on k{r) of the inte- 
gral term in Eq. (2). The mismatch may be marginally 
alleviated if one relied on the X-ray data from ROSAT 
instead of XMM-Newton that with its higher-resolution 
instruments may enhance the dumpiness effects, so bias- 
ing high the brightness (see Churazov et al. 2012). On 
the other hand, we recall from § 1 that an analogous SZ 
vs. X-ray mismatch is obtained with quite different fit- 
ting tools by Ade et al. (2012). Thus the tension turns 
out to be model-independent, and calls for a physical ex- 
planation in terms of a pressure contribution adding to 
the thermal value p. 

In this context, the SM formalism is endowed with 
an extra gear (Cavaliere et al. 2011), i.e., its ability 
to straightforwardly include in the equilibrium a non- 
thermal component 5p to yield the total pressure p + Sp. 
The result can be simply described in terms of Eq. (2), 
with p and k rescaled to 

p = p{l-\- 5p/p) , k = k{l + 5p/p) . (3) 



In the above paper we have discussed in one particu- 
lar instance (focused on turbulence in cluster outskirts) 
how the pressure 5p{r) can be physically characterized 
in terms of a normalization provided by the infall kinetic 
energy seeping through the virial shocks to drive turbu- 
lence, and of a dissipative decay scale. 

For Coma a decay scale is not needed, and a nearly 
uniform Sp/p applies to a good approximation. Thus the 
net outcome is to lower the normalization applying to the 
thermal pressure at the virial radius, to read pr (x {I + 
Sp/p)~^. Resolving the tension between the SZ vs. the 
X-ray data requires 6p/p « 15% (up to 20% for XMM- 
Newton, which may however include a 5% bias due to 
dumpiness, see above). The outcome is illustrated in 
Fig. 2 by the solid line; we remark that while the SZ 
profile from the SM has not been derived from a formal 
fit, yet it turns out to represent well the Planck data 
over their whole radial range. In particular, the thermal 
pressure derived with the SM is now lower by ~ 15%, as 
in fact sensed by the SZ effect. 

Note that a uniform Sp/p implies the density n oc 
[p/k)^/^ to be closely unaltered, while the temperature 
normalization is affected as Tr c>c prJur oc {1 -\- Sp / p)~^ ; 
this amounts to a minor recalibration in the strength of 
the virial shocks (see § 2; also Cavaliere et al. 2009; 
2011). The resulting temperature profile stays close to 
that in Fig. 1 (bottom panel) with the thermal compo- 
nent of the central entropy recalibrated to k^ ~ 470 from 
the previous value around 540 keV cm^. 

To sum up, the thermal electron pressure is related to 
the equilibrium pressure p by 
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Figure 3. Pull spectral distorsions of the CMB intensity due to the SZ effect, computed for the Coma cluster. The lower scale represents 
the quantity x = hu/kgTQMB) while the upper scale is labeled with the frequency u in GHz; note that Planck is sensitive to bands in 
the range u 100 — 857 GHz. Magenta dotted line refers to the thermal SZ effect; red solid line is the thermal effect with the relativistic 
corrections for the Coma average temperature fc^T = 8.2 keV; blue dashed line is the non-thermal SZ effect from relativistic electrons 
down to 71 = 10^ with density 10~^ cm~^, and green dot-dashed line down to 71 = 1 with density 10~^ cm~^. Note in the latter the high 



frequency tail, and the displacement of the null from the thermal value, see 
spectral index s = 3.4 is adopted. 

With Sp/p w 15 - 20%, this boils down to pe w 0.45 - 
0.42 p, definitely lower than the bound 0.5 p pointed out 
in § 1. Note that sensible variations in the ICP metal- 
licity Z « 0.4 ± 0.03 measured in Coma by Sato et al. 
(2011) would bias only by a few percents the electron 
pressure inferred from the X-ray bremsstrahlung radia- 
tion, as discussed by Churazov et al. (2012). 

4. DISCUSSION AND CONCLUSIONS 

Next we discuss the physical nature of such a non- 
thermal pressure contribution Sp to the overall equilib- 
rium. 

• Cosmic-ray protons potentially constitute attrac- 
tive contributors (e.g., Pfrommer et al. 2005), as 
their energy is longlived and can be stored within 
a cluster. However, in Coma their overall energy 
density has been bounded to be less than a few • 
10~^ of the thermal pressure by radio and 7-ray 
observations (see Ackermann et al. 2010; Bonafede 

et al. 2012). On the other hand, cosmic rays may 
play a role as injectors of secondary electrons, to be 
subsequently accelerated by turbulence and shocks 
in the ICP (see Brunetti et al. 2012). 

• Ongoing turbulence originated by recent mergers 
that drive plasma instabilities in the weakly mag- 
netized ICP constitutes an attractive contributor 
in view of its direct link to the primary energet- 
ics. Such a turbulence has been discussed by many 
authors as a source of velocity and density fluctu- 
ations (see Nagai et al. 2007; Vazza et al. 2010; 



4 for details. A powerlaw electron energy distribution with 

lapichino et al. 2011); it is widely held to accel- 
erate with moderate efficiency supra-thermal elec- 
trons in the plasma to mildly relativistic energies 
giving rise to steep distributions (see Schlickeiser 
et al. 1987; Sarazin & Kempner 2000; Blasi et 
al. 2007; Brunetti & Lazarian 2011). However, in 
Coma the density fluctuations caused by ongoing 
subsonic turbulence have been constrained by Chu- 
razov et al. (2012; see their § 5.2 and 5.3) to be less 
than 5% on scales 30 — 300 kpc. The correspond- 
ing indirect estimates of current turbulent veloci- 
ties < 450 km would fall short of providing the 
additional pressure required to relieve the SZ vs. X- 
ray tension. The actual turbulence velocities will 
be directly probed with the upcoming ASTRO-H 
mission (http : //www. astro-h. isas .jaxa.jp/). 

Relativistic electrons with Lorentz factors 7 > 10^ 
in the diffuse magnetic field J5 w a few /iC mea- 
sured in Coma emit the large-scale synchrotron ra- 
diation observed at > 30 MHz in the form of the 
classic Coma radiohalo, see Govoni et al. (2001) 
and Brunetti et al. (2012). Based on the halo 
shape discussed in the last reference, the pressures 
of the magnetic field and of the energetic electrons 
appear to be effectively coupled to that of the domi- 
nant thermal population (see discussions by Brown 
& Rudnick 2011 and Bonafede et al. 2012). The 
integrated radio power of several 10^° erg s 
plies a relativistic energy density of order 
erg cm~'^ (see Giovannini et al. 1993, with param- 
eters updated). Although the corresponding pres- 
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sure value is substantially smaller than the required 
5p w 0.15p ~ several 10^^^ erg cm^"^, relativistic 
electrons can provide interesting candidates if their 
energy distribution extends steeply toward a lower 
end 71 < 102. 

Such an extension is consistent with the radio spec- 
trum retaining a slope a « 1.2 or somewhat steeper, as 
observed down to frequencies v w 51 MHz (see Hcnning 
1989); the corresponding electron distribution is to rise 
toward low energies as ^'^'^ with slope s = 2a + 1 ~ 3.4. 
Existing data (cf. Henning 1989) also show that at lower 
frequencies the radio flux in Coma is still sustained, and 
may even feature a steeper component, as found in other 
clusters (see van Weeren ct al. 2012); LOFAR will soon 
clear the issue (see http : /www. lof ar . org/). 

The amount of non-thermal pressure implied by the 
above electron population may be estimated as Sp w 

iTLpC^ n^icii"'!} /S oc 7^"*, and refined with the full ex- 
pressions including mildly relativistic electrons as given 
by EnfiUn & Kaiser (2000, their Appendix A). Using the 
value 2 x 10^° erg s~^ of the radiohalo luminosity at 100 
MHz and the profile given by Brunetti et al. (2012), we 
compute that a non-thermal contribution 5p/p m 15% 
would indeed obtain on extending a straight electron dis- 
tribution down to 7i ^ a few. 

On the other hand, a slope sustained against the fast 
Coulomb losses (e.g., Sarazin 1999; Petrosian & East 
2008) requires such electrons cannot be drawn from the 
thermal pool, but rather to have been injected over a 
few 10^ yr by the action of mergers, or by AGNs (like 
the current sources associated with NGC 4869 and NGC 
4874), or by cosmic-ray interactions (see Brunetti et al. 
2012). These electrons are widely held to be acceler- 
ated via turbulence and low-A^ shocks, recently driven 
by mergers deep and energetic but already on the way of 
dissipating, so as to meet the constraints set by Chura- 
zov ct al. (2012) and recalled above. We have stressed at 
the end of § 2 that similar merging events over timescales 
of Gyrs are independently required for providing the top 
level kc « 500 keV cm^ of the central entropy in Coma. 

Energy distributions steep down to 71 ~ a few im- 
ply a density n ^ 10^^ cm"'^ in trans-relativistic elec- 
trons, and so provide an upper bound for gauging the 
actual I0W-7 electron population via the tail of the SZ 
effect at very high frequencies > 1 THz, and the dis- 
placement of the thermal null at 217 GHz; these spectral 
distorsions are illustrated in Fig. 3 (sec also Rcphaeli 
1995). Such features in the SZ spectrum are within 
the reach of sc^usitive instrunic;ntations like ALMA (see 
http : //www . almaobservatory . org/). 

In the range 7 < 10^ the electron distribution will 
be progressively flattened down by Coulomb interactions 
over timescales < 10^^ Gyr. But a 'silent pool' of cool- 
ing electrons with 7 ^ 10^ can be replenished and piled 
up since their lifetimes top at about 1 Gyr. With a cu- 
mulative density n ~ 10~^ cm""^ resulting from many 
mergers, these electrons can yield a non-thermal contri- 
bution 6p/p ~ 15%. Their synchrotron and relativistic 
bremsstrahlung radiations would easily escape detection 
(Sarazin 1999; Sarazin & Kempner 2000), while their col- 
lective contribution to pressure is probed just by the ther- 
mal SZ effect. Note that to sustain such a silent pool 
would require a rather high energy dissipation rate by 



mergers. 

In siimmary, the intriguing physical conditions fea- 
tured by the inner IGF in the Coma cluster include both 
a thermal and a non-thermal component, to be probed 
via three observational channels across the electromag- 
netic spectrum: the bremsstrahlung emission in X rays, 
the thermal and relativistic SZ effects in microwaves, and 
the diffiisc synchrotron radiation in the radio band. On 
comparing the first two views, we have found a model- 
independent tension between the SZ Planck data and the 
X-ray observations. In fact, a similar mismatch has been 
also found by Ade ct al. (2012) with their empirical flt- 
ting formulae, concerning not only spherically-averaged 
profile, but also 3 out of 4 angular sectors probed in de- 
tail. 

On large scales, the tension is hard to explain out in 
terms of overall asphericities given their limited impact 
in Coma, as discussed by De Filippis et al. (2005). On 
small scales, the narrow shock jumps reported in selected 
sectors by Planck (Ade et al. 2012) are diluted by line- 
of-sight projection and azimuthal averaging. On inter- 
mediate scales of order several 10^ kpc, the presence in 
the ICP of substructures (see Ade et al. 2012) and fluc- 
tuations (see Khedekar et al. 2012) may contribute to 
locally bias the X-ray pressure; however, given the con- 
straints on the density fluctuations in Coma (Churazov 
et al. 2012), we expect these effects to be limited to 
about 5%. 

We have instead proposed that the SZ vs. X-ray ten- 
sion can be resolved in terms of a physical condition, 
i.e., non-thermal support Sp/p ^ 15 — 20%i yielding a 
lower thermal electron pressure Pe ~ 0.52p/(l -|- dp/p) 
0.45 - 0.42 p, see Eq. (4). The spherical SM proceeds to 
provide a pressure profile that straightforwardly incorpo- 
rates the non-thermal contribution to pressure, and fits 
in detail the SZ shape while being consistent with the 
resolved X-ray obscrvablcs (sec Figs. 1 and 2). We stress 
that, at variance with the fitting formulae used by Ade 
et al. (2012), the SM pressure profile also features a slow 
decline for r > 0.3 i? in agreement with the Planck SZ 
data. Related features are also pleasingly apparent in 
Fig. 1 of Lapi et al. (2012), where the SM is found to 
perform considerably better than the fitting formulae by 
Arnaud et al. (2010) in reproducing the stacked SZ pro- 
file from several clusters observed with the South Pole 
Telescope. These successes support our view that the 
discrepancy between X rays and SZ effect in Coma is not 
dominated by specific asymmetries. 

Given the current constraints on cosmic-rays and on 
turbulence (pending direct measurements of the turbu- 
lent velocities), we have discussed the additional non- 
thermal pressure in terms of a mildly relativistic electron 
population with 7 in the range from a few to 10^. We 
stress such a component to constitute a natural byprod- 
uct of the intense merger activity independently required 
for yielding the very high central entropy in the thermal 
intracluster plasma of the Coma cluster. 
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